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We build an exact description of out-of-equilibrium fixed points in quantum impurity systems, 
that is able to treat time-dependent forcing. We then show that exact analytical out-of-equilibrium 
results can be obtained in interacting quantum impurity systems in their strong coupling regime, 
provided they are integrable at equilibrium and they only allow for integer charge hopping at low 
energy. For a large class of local operators, including currents (of charge, energy, ...), we show that 
out-of-equilibrium correlators can be obtained as equilibrium correlators of effective operators. We 
explicitly show that out-of-equilibrium properties are universal, i.e. they depend on the microscopic 
details through a single energy scale. We apply our approach to the Interacting Resonant Level 
model, and obtain the expansion of the universal scaling function for the charge current as a function 
of voltage, temperature, and frequency, up to order seven. 



The description of out-of-equilibrium properties of in- 
teracting quantum many-body systems is certainly a 
major challenge of modern physics. Amidst the im- 
mensely large class of such systems, quantum impurity 
systems (QIS) are remarkable in that they can be con- 
sidered as the simplest ones: they are defined as (or can 
be brought to) one-dimensional systems homogeneous 
in space (metallic quantum wires that play the role of 
baths) , except at one point where the interaction is con- 
centrated. Moreover, they can easily be forced out-of- 
equilibrium by coupling the wires to different external, 
possibly time-dependent, sources (of charge, energy,...). 

QIS model transport through microscopic (e.g. quan- 
tum dots [ll,[3) or nanoscopic (e.g. atoms or molecules0- 
0) objects connected to macroscopic electrodes, a rich 
experimental field driven amongst other things by long 
term efforts towards the miniaturization of electronics. 
Understanding transport in QIS is thus of crucial im- 
portance. Whereas linear transport - that boils down 
to equilibrium properties - is fairly well under control 
theoretically, non-linear effects are significantly harder 
to predict: solving the out-of-equilibrium theory of QIS 
unfortunately turns out to be a considerable problem, 
mixing many-body aspects (interactions make it compli- 
cated) with the intrinsically open geometry of the out- 
of-equilibrium problem. 

The one-dimensional character of QIS comes along 
with a realm of powerful methods - be they analytical 
or numerical - that gives the hope one could solve the 
out-of-equilibrium problem. A generic feature of QIS at 
equilibrium is that the impurity/wire coupling, no matter 
how small it is, has drastic consequence on the ground- 
state of the system: at temperature T = 0, properties 
are described by a strong coupling (SC) fixed point 
that physically corresponds to the hybridization of the 
wires and governs the low voltage behavior: for example, 
for systems with a Fermi liquid SC fixed point hy- 
bridization results in a linear I{V) characteristic at small 



mal at T = 0. A typical energy scale, called here the hy- 
bridization temperature Tb (or equivalently a hybridiza- 
tion length Ab = fc^y- with Up the Fermi velocity) and 
akin to a Kondo temperature, marks the crossover to the 
SC regime: as soon as E/T^ acquires a finite value, E be- 
ing any energy scale at which the system is probed, addi- 
tional many-body scattering mechanisms must be taken 
into account to describe the physics - and, incidentally, 
non-linear effects. This SC regime, defined by all energy 
scales being smaller than Tb, is the focus of this work 
(FigllJ ; it is not captured by conventional perturbative 
expansions, performed around the weak coupling (WC) 
fixed point, that encounter convergence problems when 
the largest [i^l of the ratios E/T^i becomes < 1. 



WC 



FIG. I: A sketch of the renormalization group trajectory 
of the QIS interpolating between the decoupled (WC) and 
hybridized (SC) fixed points. The energy scale Tb marks 
the crossover between the weak coupling regime and the SC 
regime (our focus, shaded area). 

Numerical techniques have made significant pro- 
gresses recently, with the development of efficient time- 
dependent algorithms. Nevertheless, real-time numeri- 
cal approaches performed on finite systems like time- 
dependent DMRG face a difficulty in extrapolating to 




large sizes or times jllj , and methods formulated in infi 



voltage, with the conductivity Go 



di 
dv \v=o 



being maxi- 



nite systems like diagrammatic Monte Carlo [1^ or time 
dependent NRG [l^l run into problems for accessing the 
stationary state in the SC regime where the equilibration 
time becomes too large. Convergence of real-time renor- 
malization group is problematic in the SC regime. 

On the analytical side, since many QIS are integrable 
(and thus exactly solvable) at equilibrium, one could 
hope for a full exact solution out-of-equilibrium. While 
this is in fact the case for free QIS (by free we mean 
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that the hamiltonian is quadratic in fermions, a very 
special subset of integrable QIS) where the Landaucr 
Buttiker formahsm aDphes jlsMlTj . and also for some in- 
teracting QIS that are amenable to free QIS by possi- 
bly complicated, non local transformations 18l42l| , this 
might well be a property of free fermion systems. In- 
deed, there are to the best of our knowledge only two 
examples of interacting QIS that could be solved exactly 
out-of-equilibrium using integrability: the Boundary Sine 
Gordon model, solved by algebraic metho ds l 22| or by 
Thermodynamical Bethe Ansatz (TB A) [2^, |24|. and the 
Interacting Resonant Level model (IRLM) at its self dual 
point using TBA 
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25| . It is likely that the integrable 
structure is not preserved when the system is coupled to 
reservoirs, except in a few exceptional cases. Moreover, 
the exact TBA approach is limited to static forcing, and 
cannot treat e.g. the experimentally relevant AC regime. 

In this Letter, we are going to use integrability in a 
weaker sense, by showing that it can provide an exact, 
universal, and systematic expansion of physical quan- 
tities in out-of-equilibrium conditions (including time- 
dependent forcing) in the SC regime, that can be pushed 
to arbitrary order in principle. Precisely, we compute 
the universal scaling functions determining the physical 
quantities (e.g. / = V f{^, ...) for the current) 

in an expansion in T~^, building a kind of Sommerfcld 
expansion in interacting systems. We implement our ap- 
proach in Fermi liquids, but it goes far beyond the stan- 
dard Fermi liquid (FL) approach [l3| that only fixes 
the first correction to the T = linear regime: indeed, 
at orders higher than T~^, new scattering processes, not 
considered in the FL approach, come into play and affect 
non-linearities in the I{V) characteristics. In contrast to 
what happens at higher dimensions, (26l.[27l| we find that 
physical quantities are analytic. 

We proceed in three steps: First, the forcing out of 
equilibrium (that can be dynamical) is exactly incorpo- 
rated in the description of the SC fixed point for arbitrary 
QIS: we build the Hershfield's operator [i^ that deter- 
mines the out-of-equilibrium density matrix. 

Second, we focus on integrable QIS and carry on a 
Keldysh expansion in the distance to the SC fixed point. 
In the case of a super Fermi liquid (to be defined below) , 
analytical properties allow for an exact expression for 
the out-of-equilibrium average value of local operators 
A{x, t) in terms of equilibrium average value of effective 
operators in a free theory Hq'^': 



(A(2:i,ti)...)„.,q = (^f(xi,ti)...)^,, 



(1) 



where the operators undergo two operations : a dressing 
'^B by all scattering processes with the impurity, and a 
gauge transformation '%.Eq that translates the out-of- 
cquilibrium forcing. We give in Eqs. (|4l5l) the explicit 
form of these super-operators. 



Third, we obtain an explicit expression for the effective 
operators ([1]) in a systematic expansion around the SC 
fixed point, that organizes as a series in integer powers of 
where E is any energy scale in the problem: voltage, 
frequency, temperature(s)... We carry out this expansion 
at order (^)~^ summing up ~ 1500 diagrams using a 
Mathematica code. In full generality, the value of the 
radius of convergence of such an expansion is an open 
question. In the case of the IRLM investigated below, 
the radius is finite when there are no interactions, even 
in the case of harmonic forcing 2^ 30[ , and also at the 
self dual point in the static limit [25|, so that it gives 
good confidence that it remains finite for any value of 
the interaction at least for harmonic forcing. It might be 
in other situations: our expansion would then be an 
asymptotic series. 

We now implement this method in the concrete ex- 
ample of the IRLM 3l|, that has recently earned the 



status of benchmark in out-of-equilibrium physics in 
QlSdil, m, m, m as one of the simplest QIS allow- 



ing for non equilibrium in the presence of interactions. 
The IRLM consists of two baths of free spinlcss electrons 
coupled via tunnel hopping to a single level, that inter- 
acts capacitively with the electrodes via a short range 
potential with strength U. After the standard steps of 
linearizing around the Fermi points and unfolding [3J], 
the two semi-infinite wires are described by two right- 
moving free fermionic fields "01(2) (a;,t) coupled at .t = 
to the impurity level (with creation operator d^): 

/oo 
dxi^^d^i; (2) 

= -ia^l{Q)d + \i.c. + U:4,li:,:{Q){dU-\)+edSd, 

A sum over a ~ 1,2 is implied in Hb , and in the fol- 
lowing we set = e = = fcs = 1. For U = Q, one 
recovers a free theory, the Resonant Level Model. In 
the interacting model f7 7^ 0, at equilibrium, standard 
manipulations (see e.g. Ref. 35[) lead to the conclu- 



sion that the tunneling term has a scaling dimension 
D{U) = i + (■ ^■■'=tan(^^-f ) ^2^ implying that it is relevant 
[D < 1) in the entire repulsive regime U > and part 
of the attractive regime. The free theory corresponds to 
D = ^. In the following we restrict our attention to the 
interesting situation _D < 1 ; in this case the tunneling 
amplitude 7 (one sets 71-1-172 = 76*^^^) flows to SC under 
the renormalization group and reaches the value 1 at an 



energy scale Tb ~ 7 



1/(1-/5). 



below this scale, the system 



enters the SC regime (Fig[T]). 

The physics in the SC regime is governed by the SC 
fixed point, that is reached at vanishing energy ^ — ^ 0. 
Physically, this fixed point corresponds to a perfect hy- 
bridization of the wires: the IRLM at the SC fixed point 
can be described in terms of new fermions 5'a(x,i), the 
SC fields, that do not feel at all the impurity: the sys- 
tem is homogeneous, with Hamiltonian Hq'^ = Ho[^i] + 
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Ho[i'2]- There is a simple relationship between the SC 
fields "ila and the original physical fields ipa- for incom- 
ing fields (i.e. in the region x < 0) their densities are the 
same, :'0lV'f,: < 0,i) =:*3,*^: {x < 0,t). For outgo- 
ing fields (region a; > 0), the relationship is still linear, 
iV-i^f,: {x > 0,t) = Bli -.^I'^d-. (x > OJ) with B(D,e) a 



4x4 matrix that can be worked out [35 1. 

(i) Out-of- equilibrium SC fixed point: linear regime. 
We first study the SC fixed point (i.e. we set Tb = oo), 
and we force the system out-of-equilibrium by consider- 
ing a bath with thermodynamical variables {fia,Ta) that 
in turn specifies the density matrix for incoming states in 
wire a. Since densities of the original and SC incoming 
fields coincide, and since SC fields ^ are right-moving 
fields of an homogeneous system, we can prove that the 
density matrix at the SC fixed point out of equilibrium 
is pI% = pI%,i ® pI%^2^ with : 



'^N.Eq.a 



— e 



(3) 



In other words, the Hershfield's operator [28| at the SC 
fixed point is / dx fia ■'^\'^a- i^)- This seemingly simple 
operator becomes a complicated, non-local object that 
mixes the wires once reexpressed in terms of the physical 
fields Ipa- Our Eq. ([3]) generalizes an expression ob- 
tained earlier for 9 = tt/2, = and Ti ^ [s^ [H. 
Besides, we generalize this approach for time-dependent 
forcing p,a{t), where Pa(t) is an effective chemical poten- 
tial close to the impurity [i^ . Notice that as far as prop- 
erties close to the impurity are concerned, one can equiv- 
alently couple system to a space- varying chemical poten- 
tial and replace (lait) — )■ p.a{t — x) [45|. Those "right- 
moving" potentials preserve the holomorphic dependence 
of the fields, and can be absorbed exactly by a gauge 
transformation ^ a(z) %~iq{z) ■ '^aiz) = e^"<^)*a(2:) 
where Sa(a;,i) = i ^ dt' iia{t') is continuated to the 
complex plane z = i(t ~ x). The gauge transformation 
maps the out-of-equilibrium theory onto an equilibrium 
one : for a generic operator built with fermions fields, 
one has = with ^^'^ ^ ^^.^q • * 

and ( ■ ) is the equilibrium average with density matrix 

The explicit form of this transformation reads: 



^N.Eq(^) = 7^e-*4*"-"('")^*"*-(-) 



(4) 



where tzJ/'^^'-'^^ = 1+J2^^^^ ^§^dzi...§^dzkX{zi)... 
X{zk) is an ordered exponential with the contours Ck 
radially ordered around z, z] > ... > l^^fe— z| >0. Our 
Eq. dl]) proves that forcing the system out-of-equilibrium 
at the SC fixed point thus amounts deformation of the 
underlying conformal field theory describing the SC fixed 
point. This results holds irrespective of the nature of SC 
fixed fixed point - including non Fermi liquid fixed points. 

Using '^N.Eq, let us now evaluate the electrical cur- 
rent operator on wire a close to the impurity, la ='■ 



^I'^a- ~ ■'>Pl'>Pa' (0^). '^N.Eq mixes descendent 

fields amongst each other, and we find that the current 
acquires a part proportional to the identity: I°^{t) ^= 

EUo^ait) (0,0+^^(i) With F = Mi-M2 the 

voltage across the impurity. One thus recovers the linear 
regime If = (/,)„.Eq = Go V{t) with Go = sin^ 9 ^. 

(ii) Around the SC fixed point. We now consider fi- 
nite values of ^ , ^ , that drive the system away from the 
SC fixed point. The inhomogeneity of the system reap- 
pears: the SC fields 'ia now undergo scattering processes 
at a; = 0, and we should consider all irrelevant processes 
allowed by symmetries. The IRLM is a Fermi liquid, so 
the lowest order process is an energy-momentum tensor, 
call it O2, with coupling oc It determines the back- 
scattered current I^s = la' — {la) V^/T^ at lowest or- 
der: this is Fermi Liquid Theory. At higher orders T~" 
(e.g. to compute the first finite temperature correction 
cx r^y^/Tg to /bs), O2 is not sufficient and indepen- 
dent, higher order processes must be considered. In gen- 
eral one faces here several difficulties, that prevent any 
practical calculations: the number of independent pro- 
cesses (and corresponding couplings) grows very rapidly 
with n ; moreover there is no way to fix the couplings 
[iij . Therefore, physical quantities depend on a whole 
series of energy scales (the couplings are dimension-full) 
that are generated by the finite distance to the SC fixed 
point. On the contrary, as we shall see now, in an inte- 
grable system, out-of-equilibrium physical quantities de- 
pend on a single energy scale Tb, and are therefore given 
by universal scaling functions. 

The IRLM is integrablc, and its infinity of conserved 
quantities 02n highlyconstraints the allowed processes. 
It has been shown [37| that a dual description of the 
Hamiltonian can be derived: H = H^'^ + H^'-' with 



Hb'^' — J2n=l rr^Yr,-! '^2nK 

arc known exactly. This dual Hamiltonian has been used 
to derive equilibrium properties [i^lsl] in the SC regime. 
The operators 0211 are the only independent allowed pro- 
cesses, and transfer integer charges across the impurity: 
we call such a system a "super Fermi liquid" . 

We now evaluate a general out-of-equilibrium iV— point 
correlators of the currents. We perform a Keldysh expan- 
sion in powers of T^^^, starting at time t = —00 at the 
out-of-equilibrium SC fixed point (Ab = 0) and adiabat- 
ically turning on AB(i) to a finite value Ab(0) = ^ > 0. 
Using a super-operator formulation [s^ , and the fact that 
in a super Fermi liquid, nested commutators of H^'^ with 



92 



^02nix = 0) where the numbers g2n 



the currents can be represented as contour integrals 47[ , 
we arrive at the following form for the super-operator 
that implements the Keldysh expansion: 



(5) 



Importantly, short-distance divergences in the expan- 
sion arc exactly canceled by the regularization (point 
splitting) inherited from integrability, even out-of- 
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equilibrium. This proves formula ([T]), an expression that 
is finite order by order. For incoming fields, we find that 
■ "^aix <0) ~ ^a(a^ < 0) as expected from causality. 
(in) Predictions. Calculating the effective cur- 
rent operator 1°^ = X)^o exploiting ([T]) 
in a perturbative way requires the computation of 
expressions '^n.eq ■ 02711 * ■■■ * ^211^* ■ '^\^b ■ where 
O -k ^{z) = § dwO{w)^{z). It becomes a purely 

algebraic problem: at order Tg"", lives in the 

space 8n+i of Kac Moody descendent fields [i^ with 
conformal weights <n+\ At large order, we resort 
to a Mathematica code (a calculator in £n+i, a space 
with intricate structure, in which basic operations are 
neither commutative nor associative). We performed 
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FIG. 2: [Color online] Differential conductance Gbs{V) at T = 
2 

0,0 — ^ in ^ units. We also plot the exact results known at 
D — ^ (free fermions) and D = j (self-dual point). Dotted 
lines are the FL results (up to Tb"^)- The inset shows the 
temperature dependence of a, the first non-linearity in — 

a{T)^+0{V^). 

calculations up to order T~'^ , and we give here the 
expression of the first terms of the backscattered current 
-^Bs = ^0*^ — {!) for a static applied bias V, atT ^0 and in 
the resonant case = (from now on one considers the 
symmetrized current I = — 



X(V^ + (2ttTY) 

2(2irT)^ (l+5X+X'^)+W'^ {X'-' + 1 

and 54 - •^■i-x' y 



3g4(2(27rr)-'X(X-15)+3V^(X^-10X+5)) 



D 
67r^ 



+ o{^) with x = 4i:i-i 

— -. The resulting back- 



scattered non- linear conductance Gbs(T^) = is shown 
in FiglH We see the excellent agreement (up to 0.7 Tb) 
with the exact expressions available in the free and 
self-dual cases: it improves on the FL results (drastically 
close to the self-dual point D = j). 

We then consider an AC bias /Xi/2(i) = ±V/2cosLut. 
In the small voltage regime, the current can be writ- 
ten in terms of the admittance Yy'\u!) — [G + iujC){uj) 
(C is the capacitance), such that I{t) = V ^{Y'^^h"^^) , 
and non-linear effects now lie in the w-dependence. We 
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FIG. 3: [Color online] (0 = ^) (a) Modulus of the admittance 
|y(a;)| for an AC bias (in /h units). The curves are dashed 
when the last two orders of our expansion differ by more than 
5%. (b) Capacitance C(aj — 0,T), normalized by the FL 
value C(0,0). (c) Amplitude of the 3a; harmonics y'"^' as a 
function of frequency (r = 0) (d) Influence of temperature on 
the efficiency of harmonics production. 



show in Fig|3] the amplitude jH'"^''!' it is maximal in 
the non-interacting case. On top of the T, w — > limit 
C(0, 0) = predicted by FL theory, we find a rich 



4DTBh 

dependence on uj and T, in particular C(a; = 0) is signifi- 
cantly affected by temperature in a way that depends on 
the interactions. 

We are also able to compute explicitely higher 
harmonics of the current. We define higher 
admittances as Fourier coefficients ll^"^(i^) = 

limv^o ^^^^f^ Jli'^/^ I{t)e-""^*dt. In the particle- 
hole symmetric case = 0, the current is odd in V and 
the first non linear contribution is an odd harmonic at 
3a; . In Figl3]we show how frequency affects the ampli- 
tude Yj^l The FL result = ^^^~q^^^T' ^ cannot 
capture the sign change of the a; dependence when 
interactions are switched on. If we break particle-hole 
symmetry with a grid potential edd^d, 2a;— harmonics are 
generated, and the lowest contribution to the admittance 
reads y^^^ = -ig4 %l°^f ^ ■ it docs not depend on the 
first process O2, and is thus missed by the FL approach. 

Conclusions. We have shown that out-of-equilibrium 
QIS at their SC fixed point can be described by a simple 
deformation of the underlying conformal field theory de- 
scribing the equilibrium SC fixed point, in the same spirit 
as finite temperature can be obtained as a mapping from 
the plane to the cylinder geometry. Then, in integrable 
QIS that are super Fermi liquids, we have shown that the 
integrable structure can be used to build a super-operator 
mapping the out-of-equilibrium interacting theory onto 
a free one at equilibrium, and allowing us to show that 
physical properties are universal. This also results in the 
possibility to obtain controlled, exact analytical results 
in the whole SC regime for the expansion of the univer- 
sal scaling function of the currents. We have applied this 
method to the IRLM, revealing some non-linear effects 
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that are not described by the usual Fermi liquid the- 
ory. Our approach furnishes a mean to investigate in a 
controlled way many other out-of-equilibrium properties 
of super Fermi liquids, like (finite frequency) noise, full 
counting statistics, or non-linear thermal transport. It 
also raises the question whether it can be generalized to 
non (super) Fermi liquid integrable QIS. 
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